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Abstract 

We study thermalization of strongly coupled gauge theory plasma in the presence of magnetic 
held using the AdS/CFT correspondence. We utilize the falling energy-shell model as a holographic 
description of gauge theory plasma undergoing thermalization, and find the effect of magnetic held 
on thermalization time in various space-time dimensions. Our results demonstrate that magnetic 
held universally hastens thermalization of strongly coupled gauge theory plasma. 
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I. INTRODUCTION 


Recently, studying the effects of magnetic field on the properties of quark-gluon plasma 
(QGP) has received a considerable amount of attention, as the QGP created in ultra- 
relativistic heavy-ion collisions may be subject to a huge coherent magnetic held produced 
by many spectator nucleons [1]. Magnetic held can probe several key dynamical and 
topological M properties of QGP related to havor symmetry of QCD, which are not easily 
accessible by other gluonic observables, providing important complementary characteriza¬ 
tion of the created QGP. In this paper, we study the effect of external magnetic held on the 
thermalization of QGP at strong coupling regime using the AdS/CFT correspondence [7]. 
Magnetic held can potentially be important in the thermalization of QGP in heavy-ion col¬ 
lisions, since the thermalization occurs at an early stage of heavy-ion collisions when the 
magnetic held is strong, before it dies out with time. 

In our study, the QGP that undergoes thermalization is modeled by falling of a thin spatial 
mass shell to the bottom of AdS space, forming a black-hole at the end of thermalization |[8j. 
See Refs. [91 IT4] for other approaches. We treat the magnetic held as external, and use 
known solutions of AdS geometries with magnetic held at zero and finite temperatures. In 
a thin-shcll approximation, we join two static solutions, one with zero temperature and the 
other with finite temperature, across the falling shell via the Israel junction condition [15j. 
At each time, the location of the shell in the energy coordinate (holographic coordinate) 
divides the AdS space into two regions: one with the geometry of zero temperature that is 
not yet thermalized, and the other with Unite temperature that is thermalized. As the shell 
falls down towards infrared, eventually forming a black-hole, the AdS space becomes filled 
with the geometry with finite temperature, representing dynamical thermalization. The 
proper time (or Eddington-Finkelstein time) by which the shell forms a black-hole can be a 
reasonable definition of thermalization time in the model [[TO] . 


II. ADS GEOMETRIES WITH ZERO AND FINITE TEMPERATURES 

For joining of two static solutions across a thin falling shell to work, each static solution 
one uses for the two different regions that the shell divides the space-time into, must be 
isotropic and homogeneous: this requirement is seen in the Israel junction condition in the 
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subsequent analysis. The underlying reason for this requirement can be understood by the 
Einstein-Maxwell equations with a source (the shell) viewed as an initial value problem. The 
shell (which is assumed to be neutral) starting from rest at initial time and moving along 
its trajectory would normally source metric perturbations inside its future light cone while 
it falls down by its own gravity. There is no a priori reason to expect that the resulting 
geometry in the future will simply be given by joining of two static geometries across the 
falling shell: one instead expects gravitational waves emanating from the shell. In the 
presence of isotropy and homogeneity however, a powerful uniqueness theorem of Einstein- 
Maxwell theory dictates that a (neutral) homogeneous and isotropic solution in a connected 
region with no sources is completely fixed by its conserved energy density, and must take 
a form of static black-hole with that conserved energy density. The two regions bounded 
by the shell have constant energy densities differing by the energy density of the shell, and 
since these energies are conserved, the geometries in each region are fixed by the uniqueness 
theorem to be those static geometries with conserved energy densities. This is the physical 
reason why the falling shell ansatz works: simply put, no gravitational radiation is possible 
in isotropic and homogeneous collapse m ■ 

We look for AdS geometries with magnetic field which possess isotropy and homogeneity. 
In dimension D = 5 (corresponding to 4-dimensional QGP) a single magnetic field nec¬ 
essarily breaks isotropy. To overcome this difficulty, we consider Af = 4 super-Yang-Mills 
gauge theory with global SO(6)r R-symmetry that allows three orthogonal magnetic fields 
from each f/(l) 3 C SO(6)r of equal magnitude. Although our model for the magnetic field 
from R-symmetry of Af = 4 super Yang-Mills theory has differences from QCD, such as the 
charge content of the matter fields, we expect that the universal feature we observe in this 
model could indicate the similar trend in real QCD in strong coupling regime. 

The corresponding theory in AdS is the gauged U( l) 3 supergravity which is a particular 
Einstein-Maxwell-Scalar theory. It admits an exact solution with three orthogonal magnetic 
fields of equal magnitude that ensures isotropy and homogeneity of the energy-momentum 
tensor (TBj. The action is given by 

(1 6ttG 5 )£ = (R - V) - \ J2 W') 2 -jE ' Y «~ 2 ( F “) 2 + r7 L=e'“" > " A . F^F%A \. (1) 

Z /-i 4 9 5 
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where F^ u = <9 M A“ — d u A c f \ and 




a= 1 


X x = e"76 01 “^ 02 , X 2 =e~^ 1+ ^\ X 3 = eVS 




( 2 ) 


We will set L — 1 in the following]^] The field equations derived from the Lagrangian ([!]) 
admit an exact magnetically charged AdSs black hole solution |18j 


2 dz 2 

ds 5 — 77 \ 9 

f(z)z 2 


f(z ) 2 (dx) 

J -±Adt 2 + W 


F“ = e aij B, 0 7 = O, 


(3) 


where a = 1,2,3 labels three 17(1) R-symmetries and i , j = 1,2,3 are spatial indices. The 
function f(z) is 

f(z) = 1 — mz 4 + \B 2 z i log (mz 4 ) , (4) 
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with zh = m~4, being the location of the black hole horizon solving f(zjj) = 0. The 
parameters (m,B) are related to the temperature T by 

f'(z H ) 8 m - B 2 


T = 


(5) 


47T 87T77T 3 / 4 

Note that at m — \B 2 , the temperature of the black hole (|5j) goes to zero, and hence the 
extremal zero temperature solution in the presence of magnetic fields is given by 


fo(z) = 1 - tfiV + \b 2 z‘ log (IfiV). 


( 6 ) 


We map this to the field theory vacuum in the presence of magnetic field. We must have 
m > |R 2 for thermodynamic stability. 

In 4-dimensional AdS space (corresponding to a field theory in 3-dimensions), one can 
realize isotropy and homogeneity with a single magnetic field F l2 = B. The exact black hole 
solution with magnetic field in the Einstein-Maxwell theory is known 

dz 2 


dsl = 


f(z)z 2 z 2 


l^) dt 2 + l Xl 


where 


2 4 


f(z) = 1 - mz 6 + B 2 z 


(7) 


( 8 ) 


1 More precisely, L 4 = ^ M AT C Z 4 but the string scale l s which serves as a UV cutoff cancels in any final field 
theory observables, and for convenience, we set L = 1. In this choice, we have G 5 = ir/(2N 2 ). 
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The location of the black hole horizon is given by f{zn) = 0, and the temperature T is 

T = f( z H) = 3 - B 2 z a h 
A'k Attz^ 

When m = m 0 = -^B 3 / 2 , the temperature of the black hole solution |9j) becomes zero, and 
hence the extremal zero temperature solution is given by the blackening factor, 

/o(.~) = l-^B 3/ y + 5V. (10) 

The metric ([t]) with m — 0, that is f(z) = 1 + B 2 z 4 , is a solution of the Einstein-Maxwell 
equation without black-hole horizon. The reason why it can not be the solution for zero 
temperature is its violation of causality: the speed of light in the bulk AdS with respect to 
the field theory coordinates (£, x) at position z is c(z) = f(z), which has to be less than 1 to 
respect causality of the field theory [20]. This means that this geometry should be excluded 
in a meaningful AdS/CFT correspondence. 


III. HOLOGRAPHIC THERMALIZATION WITH MAGNETIC FIELD 


The thin shell initially starting from rest at a position z t = 1/tiQ s collapses from the 
UV region of small z to the IR region of large z under its own gravity, eventually passing 
through its black-hole horizon by which we have thermalization. The geometry is constructed 
by joining a black hole solution with finite temperature Q above the shell in the UV region 
with the zero temperature solution ([6]) below the shell, across the trajectory of the shell in 
(t, z) coordinates that is determined by Israel junction conditions. 

The metric induced on the 4-dimensional world-volume £ of the shell can be written in 
a conformal form 


ds% — 


-dr 2 + (dx) 


, , \\2 . (ii) 

Mv) 

where z(t) is the position of the shell in z coordinate at a conformal time r. Continuity 
of the metric across the shell requires identifying x on £ with x in the background. The 
trajectory of the shell with respect to (tu, z) coordinates in the upper (UV) region of space- 
time parameterized by the conformal time r, that is (tu(r), z(r)), determines the induced 
metric on £. Comparing time component of that with 0 gives 

i 2 (r) 


/ ( Z ( T )) tu( T ) ~ 


f 0(u» 


= 1 


( 12 ) 
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where • = jk This relates tu and r, given a trajectory z(t). Similarly, the same trajectory 
with respect to the IR coordinates (£ l ( t ), 2 :( t )) should satisfy the condition 

A « t M ( t) - = 1 , (13) 

that gives a relation between and r once the trajectory z(r) is found. Finally the Israel 
junction condition is 

[Kij - 7 i3 K] = -8 t rGsSij , (14) 


where [A] = A^ — Ajj, Sij is the energy-momentum tensor on the shell, and 7 ^ is the induced 
metric on the shell with respect to the shell coordinate £* = (r, x). The K^J L are extrinsic 
curvatures evaluated on the shell from the upper and lower regions respectively, 


K ‘> = ww v ‘ n ' = ~ n ° 


d 2 x° 


d?d£ 


j + r 


Pi 


dx 3 dx 1 

ww 


(15) 


with the unit normal vectors n^/ R to the surface E pointing to the direction of increasing z. 
They are given by 


riu = 

n L = 


zz ^ + (*/(*» | ■ 


mj at 


zz 'l+(*/.woA 


(16) 


Mz))dt dz 

where all quantities are evaluated on the shell. The non-vanishing components of K^j L are 


K u = Ju_ f f(f' + 2z) _ f_ 
TT z \ 2 (/ + i 2 ) z 


K u = 
13 


i,j = 1,2,3, 


j^L = (^0 + 2 ^) /o 


K L = 
13 


z V 2 (/o + i 2 ) 2 

hfo 


.2 


5*7, i,j = 1,2,3. 


where ' = jz- Assuming energy-momentum tensor on the shell of a conformal form, 

Sij = 4 p(z)uiUj + 7*70(2:), u* = 0, 0, 0^ , 

with the pressure p(z) to be determined, the junction condition becomes 


(17) 


(18) 


fot-L—Hu 

. zfo (f + 1 ) . 2 / (f + 2 ) 

fl (/0 + 2 2 ) “ tv ~ iJ +¥)~ 


87 tG 5 p(z) , 


4 • 87rG 5 p(^). 


(19) 
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B (fm 2 ) 

z H (fm) 

m (fm 4 ) 

C (fm" 4 ) 

0 

0.209 

525.5 

263.7 

2.06 

0.209 

527.6 

262.6 

5.16 

0.208 

538.7 

260.2 

25.81 

0.187 

810.6 

272.1 

51.61 

0.162 

1469.1 

331.0 


TABLE I. (AdSs) Parameters of our numerical solutions for RHIC with a late-time temperature 
T = 300 MeV and several exemplar values of B = 0; 0.08 GeV 2 ; 0.2 GeV 2 ; 1 GeV 2 ; 2 GeV 2 . 


Removing p(z) from the above equations and using 


tL — 


\//o + z‘ 
fo 


, tjj — 


\/T+J 

f 


from (12) and (13), the resulting equation for i is integrable to give 


z = 


Cz A f 0 (z) - f(z) 

2 2 Cz 4 


fo(z) , 


( 20 ) 


( 21 ) 


with a constant of motion C > 0. This reproduces the one in Ref. [16]. We choose to express 
the falling trajectory in terms of the boundary time tjj which can be identified with the held 


theory (QCD) time on the boundary. Using the relation (12), the solution (21) translates to 


dz 

dtjj 



Cz 4 , fo{z)-f{z) \ 2 
2 ” r 2 Cz 4 J 


fo(z) 


(Cz 1 , fo(z)-f(z) \ 2 
^ 2 “ r 2Cz 4 J 


- (AW - /W) 


( 22 ) 


which we solve numerically. More precisely, the thermalization time is defined as the 
Eddington-Finkelstcin time when the mass shell passes through its black-hole horizon [16] . 

Following [16], we set our initial condition of the falling mass shell in terms of the satu¬ 
ration scale Q s , which governs the initial gluon distribution, as 


2 (tu = 0) = Zi = , z (' t v = 0) = 0 . (23) 

We measure z in units of fin. For RHIC, we take Q s = 0.87 GeV = 4.42 fm _1 , and for LHC 
we have Q s = 1.23 GeV = 6.24 fm _1 . In Table 0 we show parameters of our numerical 
solutions after fixing the final thermalization temperature to be T = 300 MeV for RHIC for 
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(a) (b) 

FIG. 1. (AdSs) Thermalization history of falling mass shell for RHIC (left) and LHC (right). The 
late-time temperature is fixed to be T = 300 (400) MeV for RHIC (LHC), and the magnetic fields 
are B = 0 (0) (solid blue); 0.08(0.3) (orange); 0.2(0.52) (green); 1(1.32) (red); 2(2.64) (violet) 
GeV 2 for RHIC (LHC). Thermalization time is when the curve reaches its plateau at the horizon. 


several exemplar values of magnetic field. In Figure [lj we show the time history of falling 
mass shell in the field theory (QCD) time tjj for a few exemplar values of magnetic field B 
with a fixed final temperature. The plots clearly indicate that the presence of magnetic field 
speeds up the thermalization of the plasma: the stronger the magnetic field, the shorter the 
thermalization time. More precisely, the thermalization time is defined as the Eddington- 
Finkelstcin time when the mass shell passes through its black-hole horizon [IB]. However, it 
is qualitatively similar to the time in Schwarz coordinate tu we show when the mass shell 
falls close to the horizon. 

Instead of fixing final temperature, we also study the case where the energy density 
measured from zero temperature but finite B state is fixed while we vary magnetic field, 
that is we fix Ae = e(T,B) — e(T = 0 ,B), which can be interpreted as the energy density 
thrown by colliding nuclei into the background magnetic field. Explicitly, we have 


N 2 

Ae = ^(^m+^B 2 


iog 1 L) 


_ , „o . , _ , , , (24) 

2 16 1 ° 1 ' ' ' ’ v ; 

which determines the parameter m in the geometry, given a fixed Ae and varying B. In 

Figure [2j we show the resulting time trajectories of energy shell with Ae chosen to be the 

energy density of T — 300 MeV, B = 0 state. Our observation of faster thermalization with 

magnetic field seems robust. 

To examine whether our conclusion depends on the number of dimensions the field theory 








FIG. 2. (AdSs) Thermalization history of falling mass shell for fixed energy density, and varying 
magnetic field B = 0 (solid blue); 0.08 (orange); 0.2 (green); 1 (red); 2 (violet) GeV 2 for RHIC. 

resides in, we study the thermalization of plasma in magnetic field in one less dimension. 
In AdSi (corresponding to 3-dimensional field theory), the analysis is the same with Q in 
the place of ([I]), and ( 10 ) in the place of Q, but with the energy-momentum tensor on the 
shell taking a 3-dimensional conformal form, 


Sij = ip(z)UiUj + 'Jijp(z) , U{ — ( j, 0, (' ) . 


(25) 


We arrive at 


dz . , 

dT n = /(z) 


\ 


(Cz £ _|_ fp(z)-f(z) 


V 2 


2Cz 3 


fo(z) 


(ai + M|Ww) 2 _(/ 0 ( 2 )-/( z) )’ 


(26) 


which can be solved numerically given the constant C which, as before, should be determined 
from initial conditions. In Figure [3j we show the time history of falling mass shell trajectory 
in field theory time tjj for a few exemplar values of magnetic field B for 3-dimensional gauge 
theory with a fixed final temperature. Again, the plots clearly demonstrate that the presence 
of magnetic field hastens the thermalization. 


IV. SUMMARY AND DISCUSSION 

In the framework of AdS/CFT correspondence, we have studied the thermalization of 
strongly coupled gauge theory plasma in the presence of magnetic field, utilizing simplified 
picture of thermalization as falling of a thin homogeneous energy-shell towards the black- 
hole horizon. Our results in various dimensions have revealed that magnetic field universally 
hastens thermalization in strong coupling regime. At weak coupling, a strong magnetic field 
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FIG. 3. (ACIS 4 ) Thermalization history of falling mass shell in AdS 4 for a late-time temperature 
of T = 300, and the magnetic fields are B = 0 (solid blue); 0.08 (orange); 0.2 (green) GeV 2 . 

causes the dimensional reduction of the system into 1-dimensional one with lowest Landau 
levels, and one may study the effects of magnetic field to thermalization at weak coupling 
in this context. It would be interesting to see how weak coupling result compares with our 
conclusion in this work at strong coupling. 

Note added: while our work was near its final stage, there appeared (2T| which addresses a 
similar question in a different set-up. Our numerical results are in qualitative agreement with 
theirs: for the values of magnetic field relevant at RHIC and LHC, the effect of magnetic 
held to thermalization time is less significant. Our work is about the universal trend of 
strong magnetic held hastening plasma thermalization. 
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